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Abstract 

We study the bulk viscosity of strange quark matter(SQM) in den- 
sity dependent quark mass model(DDQM) under the background of self- 
consistent thermodynamics. The correct formulae, with which the vis- 
cosity can be evaluated, are derived. We also find that the viscosity in 
DDQM can be higher by 2 to 3 orders of magnitude than MIT bag model. 
We calculate the damping time scales due to viscosity coupled to r-modes. 
The numerical results show the time scale can't be shorter than 10 _1 s. 
PACS numbers: 97.60.Jd, 12.38.Mh, 97.60.Gb 

1 INTRODUCTION 

Since Witten conjectured^ that the strange quark matter(SQM) composed 
of comparable number of u, d and s quarks might be absolutely stable or 
metastable phase of nuclear matter, many theoretical and observational efforts 
have been made on the investigation of its properties and potential astrophysical 
significance |21 El EI El El El- Because the difficulty of quantum chromodynam- 
ics (the lattice approach) in the non-perturbation domain , ones used to adopt 
phenomenological models. One of the famous models is the MIT bag model 
with which many authors study the equation of state(EOS) of SQM and the 
mass-radius relation of the assumed strange stars|Hl El 0. Another alterna- 
tive model is density-dependent quark mass modcl(DDQM). Originally, some 
pioneer investigations^ El OH HH E3 were carried out but the wrong thermo- 
dynamics treatment is applied. Latterly, some works dealt with the thermody- 
namics with density-dependent particle masses self-consistently with which the 
structural properties of strange stars have been re-discussed |13l ITi] . 

It had also been recognized the effects of the dynamics of SQM on strange 
star's rotations by the time that the EOS of SQM was gradually uncovered. 
Wang and Lu find that the non-leptonic weak reaction dominates bulk viscosity 
of SQM 15 . Further investigations of the bulk viscosity also have been carried 
out by many other authors in MIT bag model and the fruitful results have been 
obtained [El E3IIE1- A recent study incorporate the effective quark masses 
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in quasiparticle description into the bulk viscosity^ an d have a significant 
result for relevant astrophysics j50j. In parallel, a previous work made the study 
of the bulk viscosity in mass-density-dependent model but the self-consistent 
thermodynamics are overlooked |21). Here we will concentrate the problem in 
the background of the correct thermodynamics and consider the effect of variable 
quark masses on the bulk viscosity of SQM. As seen latter in this paper, the 
effect have twofold contributions. The first, the modification of EOS influences 
the viscosity. The second, the density-dependent masses directly contribute to 
the dynamical quantity. 

We organize this paper as such. In Sec. II, the self-consistent thermody- 
namics with density-dependent particle masses are recalled and then further 
the thermodynamical relations are discussed. In Sec. Ill, the bulk viscosity of 
SQM with density-dependent particle masses is derived, which arises from the 
non-leptonic weak interaction. In Sec. IV, the numerical results are given and a 
application to the recent astrophysical problem is discussed. In Sec. V, a brief 
summary is made. 



2 THERMODYNAMICS WITH DENSITY 
-DEPENDENT PARTICLE MASSES 

The self-consistent thermodynamics in mass-density-dependent model have been 
discussed in many works [T2l 1131 114| . However Peng et al^lj pointed that the 
thermodynamics treatments given bv[51 ll(JllllM12| have serious problem. They 
re-derived the pressure and energy expressions from the general ensemble theory 
and found an extra term can add to the expression of pressure but not to the 
energy formula. Here we don't repeat the derivation but simply write the results. 
When particle masses depend on baryon number density or the volume of the 
system, the pressure and energy have 



p= _ d(VQ(V)) 



dV 



(1) 



Although the thermodynamics potential density depends on the volume, we still 
obtain the thermodynamic relation 

E = VQ(V) + /j,N + TS, (3) 

where (3 is the reverse temperature, V is the system volume, is thermodynam- 
ics potential density, fi is the chemical potential and N is the particle number. 
According to Eq. (3) , the basic thermodynamic equation is written as at invari- 
able temperature. 

dfl 

dE = d(Vfi(V)) + d(/xA) = (0 + V—)dV + fidN, (4) 
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Therefore, the pressure have an extra term 



Equally, it is the formula ,P = — f2 + n&J^-, given in ^3]. Those previous works 
have shown that the thermodynamics are self consistent, where an extra term 
is appended to the pressure but not to the energy. In the framework of self 
consistent thermodynamics, the EOS of SQM and the structural properties of 
strange stars have been investigated. At zero temperature, the thermodynamic 
potential density is 



^ 8tt 2 



M - m^(2tf - 5m?) + 3mf In » + ^ mf 



(6) 



Accordingly, the energy and pressure can be obtained from Eqs (3) and (5), and 
the number density is 

ni = l -{£-^r 2 - (7) 

It has been found that the density-dependent quark masses should take the 
form|13| 

rrii = rriio H jyj, i=u,d,s, (8) 

where D is a parameter to be determined by stability arguments. For given 
values rib, 

n u + n d + n s 
nb = jj , (9) 

we very easily determine the quark chemical potentials /ij to evaluate the ther- 
modynamic quantities of the system by the beta equilibrium and the charge 
neutrality conditions 

fl d = fl s = U, U s = U u + fl e , (10) 

2 11 

^n u - -n d - -n s - n e = 0. (11) 

In the following, we will consider the dynamics of SQM( bulk viscosity) near 
the equilibrium state, which is governed by the EOS calculated with the self 
consistent thermodynamics. 



3 DERIVATION OF THE BULK VISCOSITY 

In a previou study [23 > the formulae in MIT model, which is presented by Wang 
and Lu^S] and developed by other investigators El EH- was simply in- 
troduced to evaluated the bulk viscosity of SQM in density dependent quark 
mass model. The investigation used an incorrect thermodynamical results with- 
out regard to the consequences due to the dependence of the quark masses on 
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baryon number density. In this section, we derive the bulk viscosity of SQM in 
the self consistent thermodynamics. Since the quark masses are the functions 
of the baryon number density represented in(6), we for convenience assume 
the baryon number density instead of the volume as had adopted in previous 
works (T51 IT1)1 IT7I I2T] oscillates periodically, 

2nt 

n b = n bQ + An b sin( ) = n b0 + Sn b , (12) 

r 

where n b Q is the equilibrium baryon number density and An b is the amplitude 
of perturbation, r is the period. The presence of viscosity results in the energy 
dissipation of pulsations in bulk matter. On one hand, using the standard def- 
inition of the bulk viscosity of matter, energy dissipation rate per unit volume 
average over the oscillation period r can be written as 

E kiD = -- / dt(divv) 2 . (13) 
T Jo 

where v is velocity of the fluid. On the other hand, the hydrodynamic matter 
deviates from the equilibrium accompanied by the time variations of the local 
pressure, P(t). The dissipation of the energy of hydrodynamic motion due to 
irreversiblity of periodic compression-decompression process can be expressed 
as 

E diss = -4r / dtP(t)V. (14) 

TV JO 

According to the continuity equation h b 4- ribodivv and the specific volume V = 
bearing in mind that Euin = — -Ediss, from Eqs(13)and (14) we have the bulk 
viscosity 

It is proven that the bulk viscosity of SQM is mainly produced by the nonlep- 
tonic interaction as follows 

u + d < — > s + u (16) 

thereby at quasi-equilibrium the pressure can be regarded as a function of the 
number density n^, where i takes u,d,s. From the baryon number density(9), 
it is suitable to assume that P = P{n b , rid, fi s )- Thus the pressure can expand 
in small deviations 

dP dP dP 

P = Po + (^-)5n b + (^-)5n d + (?—)Sn s , (17) 
on b oiid on s 

where all derivatives are taken at equilibrium and the changes in d- and s-quark 
number density will be composed of two parts 

5V 

5rii = Srii-nt—, (18) 
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where Shi denotes the number density at a given volume for s- and d-quark. 
From the reaction (16), it evidently satisfies the relation below 



5h d = -5h s = [ ^j^-dt. 

Jo dt 

Due to the baryon number conservation, we also have 



5n b 
n b 



SV_ 



(19) 



(20) 



Finally, Eq(17) becomes 



P = Po 



n b 



dP 
dn b 



dP 

dn d 



dP 

' dn s 



Sn b 


' dP 


dP~ 




n b 


dn d 


dn s _ 


I 



dh d , 

-^d, (21) 



Obviously, only the latest term contributes to the integral in Eq(15). In the 
formula above the net rate per unit volume can be transform as from the rate 
per mass given bv|21) 



= JL G% sin 2 9 cos 2 9JT 2 [ 1 



6n 2 



Att 2 T 2 



5fi, 



(22) 



where J is given in |22| and Sfi = fi s — )i d . Analogously to the pressure, the 
chemical potential difference can be expanded, giving 



an b on d an s 



(23) 



To obtain the reaction rate and then calculate the bulk viscosity, the acquire of 
the derivatives 4^- and 2^ is necessary with the equation of state in section 
above. We respectively regard the thermodynamic potential and the chemical 
potential as the functional form below 



O = rL(fij(m,nb(ni)),n b (ni)), fi, = ^i{n u n b (nj)) 
where i,j take u,d,s. From (5) and (7), we immediately obtain 



dP 

drii 



drii 



I't 



- 3n b ) 
dn b 



dnj' 



drrij 



(24) 

(25) 
(26) 



3/ijn, dn b ^ dn b 

Thus, the part in(21) contributed to the integral (15), SP, and the chemical 
potential difference at quasi-equilibrium, Sfj, can be expressed as 



SP = 



n d - 3n b fi 2 - m 2 d 



n d 



3(t 



n s — 3n b fir — mi 
n s 3/i 



(i 



5 



5jJL 



' rib(m s — md) drrid m s — m d \ Srif, ( [ir — m 



dn b 



3/i 



rib 



3rin d 



Substituting Eq(27) into Eq(15), we finally get the viscosity 



c 



1 n b o (rid — 3n & jj? 



3nb /j, 2 — m 2 s 



dt [ ^ 
o Jo at 



7r Anb \ nd 3/i n s 3/i 

For the sake of convenience, Eq(28) can be equivalently rewritten as 



dS/i 



n b (m s - m d ) dm d m s - m d \ Sn b ( /i 



li 



dn h 



3// 



rib 



SfiTld 



3/m s 



r* dh rt 



dt 

(28) 



(29) 



m 2 \ dh d 
~df' 

(30) 



4 RESULTS AND DISCUSSION 

Eqs (22) and (30) must be solved numerically to give the r&te^f- and the 
numerical result of Eq(29) can be obtained. 

In figur^we have plotted viscosity vs. relative perturbation respec- 
tively for (i)Z>3 = 156MeV, m s0 = 80MeV at n b = 0.4fm~ 3 and r = 10" 3 (solid 
curves) and (ii)Z?2 = 156MeV, m s o = 140MeV at rib = 1.36fm -3 and r = 
10 -3 (dashed curves). 




Figure 1: Bulk viscosity as a function of relative perturbation ^p*-. Solid 

curves for n b — 0.4fm -3 , = 156MeV, m s0 = 80MeV and dashed curves 
for n b = 1.36fn:r 3 , £>3 = 156MeV, m s0 = 140MeV. The both for r = 10~ 3 and 
temperature 10~ 5 , 10~ 4 , 10" 3 , 10°, 10~ 2 , lO^MeV from bottom to top. 
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In figure[21we gave the results from our formulae (solid curves) and those in 
ref.|21|f dashed curves) for the same parameters rib = 0.4fm~ 3 , D 2 = 156.0Mev, m S Q 
80Mev. Evidently our viscosity can be higher. 




-6 -4 

Log (An In ) 



Figure 2: A comparison of the our bulk viscosity(solid curves) with that in 
ref [21] (dashed curves)for parameters n& = 0.4fm~ 3 ,£>2 = 156.0Mev, m s o = 
80Mev andr = 10~ 3 . The corresponding temperatures of the curves as figure ^ 



In figureJjUwe show the results in MIT bag model from Zheng et al^Hl for 
m s o = 80MeV, r = 1CT 3 . The solid curves are because of the consideration of 
screening mass effect and the dashed curves are regardless of the coupling among 
quarks in MIT bag model. Comparing figurc[21 with solid curves in figure^ we 
notice that our results are higher by 2 to 3 order of magnitude than those of 
bag model. We also find that the solid curves in figure|21is also higher than the 
dashed curves. In a word, the facts testify to the viscosity increase of SQM due 
to medium effect. 

Recently it has been realized that rotating relativistic stars are generically 
unstable against the r(otational)-mode instabilitv|23| The viscous dissipation of 
dense matter has relevance to the onset of r-mode instability, gravitational radi- 
ation detection and the evolution of pulsars [24l 1251 1261 127j . A main motivation 
for the studies of viscosity is to calculate the damping time for compact star 
vibrations. For a star of constant density, the r-mode energy is estimated as 



E : 



2m(2m + 3) 



(m+ l) 3 (2m+ l)!^ 5 ^ 2 



The mode dissipation energy can be computed in term of 

HP f 

— = / (5<rSa*d 3 x. 



(31) 



(32) 
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Figure 3: Bulk viscosity in MIT bag for m s =80MeV , /Lt d =470MeV, r = l(T 3 s. 
Solid curves for interacting quark matter, dashed curves for ideal quark gas and 
the corresponding temperatures of the curves as figure ^ 



The damping time scale is thus 

t d = -|f = 1.31 x lOV^C 1 - (33) 

dt 

Figure 4-7show the damping time as functions of temperature for some typi- 
cal spinning strange star with periods T = 1.5ms, 3ms, 15ms and Keplian limit. 
In comparison with harmonic oscillations, the damping time scale of r-mode 
with period r = |T = 10~ 3 s(figur^J is only short to about Is, but that of 
harmonic oscillations can be short to 10 _3 s(see figure 4 in Ref 17). This implies 
that the simple radial pulsations can be damped more efficiently by the bulk 
viscosity. For a Keplian rotation star, the minimum time scale is about 10 _1 s. 



5 CONCLUSION 

Considering the self-consistent thermodynamics in density dependent quark 
mass model, we correct the error in|2J and re-derive the bulk viscosity of SQM. 
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The viscosity in DDQM is larger than MIT bag model. Consequences is that 
the r-modes of strange star will be damped more efficiently in DDQM implying 
the r-mode instability window of strange stars should be modified. 
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Figure 4: The time scales as function of temperature for a strange star 
with M = 1AMq,R = 10km, T = 1.5ms for relative perturbation 
10~ 2 , KT\ 10°, 10~ 3 , 10" 4 , 10" 5 , 10" 6 , 10" 7 from bottom to top. 
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Figure 7: As Fig 4, but assuming a Keplian rotation star. 
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